Theory of dendritic growth in the presence of lattice strain. 
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Elastic effects due to lattice strain modify the local equilibrium condition at the solid-solid in- 
terface compared to the classical dendritic growth. Both, the thermal and the elastic fields are 
eliminated by the Green's function techniques and a closed nonlinear integro-differential equation 
for the evolution of the interface is derived. In the case of pure dilatation, the elastic effects lead 
only to a trivial shift of the transition temperature while in the case of shear transitions, dendritic 
patterns are found even for isotropic surface energy. 
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Solvability theory has been very successful in predict- 
ing certain properties of dendritic growth and a num- 
ber of related phenomena (see, for example, [1, 2]). The 
solution of the two-dimensional steady state growth is 
described by a needle crystal, which is assumed to be 
close in shape to the parabolic Ivantsov solution [3]. If 
anisotropic capillary effects are included, a single dynam- 
ically stable solution is found for any external growth 
condition. This theory has also been extended to the 
three-dimensional case [4, 5]. The capillarity is a singu- 
lar perturbation and the anisotropy of the surface energy 
is a prerequisite for the existence of the solution. In the 
case of isotropic surface energy, dendritic solution does 
not exist and instead the so-called doublon structure is 
the solution of the problem [6, 7]. 

Usually, the structural phase transitions in solids are 
accompanied by small lattice distortions leading to elas- 
tic deformations (for a review, see [8, 9]). One of 
the well-known consequences is a thermodynamic elastic 
hysteresis-i.e., the splitting of the phase equilibrium point 
into two points: the point of the direct and the inverse 
transition. This is mainly due to the coherency at the in- 
terface boundary, meaning that the lattice layers remain 
continuous through the boundary. Correspondingly, the 
hysteresis disappears without interface coherence [10, 11]. 
However, the systematic investigation of the growth ki- 
netics of such phase transitions is much less developed. 
Recently, the pattern formation processes controlled by 
interface kinetics have been considered in [12] and the 
growth of spherical inclusions under elastic and thermal 
influence were investigated in [13] by means of the phase- 
field model. 

In this Letter, we discuss the influence of elastic strain 
on dendritic growth in solids controlled by heat diffusion. 
Significant progress in the description of dendritic growth 
was made by the elimination of the thermal field using 
the Green's function technique. This allows to obtain 
a closed equation for the interface evolution (see for ex- 
ample [14]). The crucial point of the present analysis is 
that the elastic field can also be eliminated by the cor- 
responding Green's function technique. By these means, 
we derive here, as in the classical dendritic growth the- 



ory, a single integro-differential equation for the shape 
of the interface which takes into account elastic effects. 
Then, we consider two simple examples of dilatational 
and shear transformations. We show that in the case of 
a pure dilatation, the elastic effects lead only to a triv- 
ial shift of the transition temperature. However, for the 
case of shear transitions, we find dendritic patterns even 
without the anisotropy of the surface energy required by 
the classical dendritic growth theory. In this sense, the 
clastic effects serve as a new selection mechanism. 

Thermodynamics of the model. Let us consider the 
growth of a new /3-phase inside of an unbounded mother 
a-phasc. We denote the characteristic lattice strain (also 
known as the stress- free strain tensor), associated with 
the phase transition, by e° fe . The free energy density in 
the initial ev-phase is 



F a = F°(T) 
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where F®(T) is the free energy density without clastic 
effects, which depends only on the temperature T, 
are the components of the strain tensor, and A and (i are 
the elastic moduli of isotropic linear elasticity. The free 
energy density of the growing /3-phase is given by: 
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Here, we neglect the difference between the elastic co- 
efficients in the two phases. We also assume that the 
elastic effects are small i.e., e° fe <C 1. Since in our de- 
scription the reference state for both phases is the un- 
deformed initial phase, the coherency condition at the 
interface reads = uf \ where Ui is the displacement 
vector. Mechanical equilibrium at the interface demands 

0nn = <Tnn and a^nr = On) , Ons = 0n*'. Here indices 
n and (r, s) denote the normal and tangential directions 
with respect to the interface; the stress tensor is defined 
as: 
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The condition of phase equilibrium requires the conti- 
nuity of a new potential 



F F 0~ nn 6 nn ^^riT^r. 
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across the flat interface [15], which takes into account the 
coherency constraint. In the general case of curved inter- 
faces, the surface energy 7 has also to be incorporated, 
and the phase equilibrium condition for each interface 
point in the case of isotropic surface energy reads 



F a - F p - 7/c = 0, 



(3) 



where n is the local curvature of the interface. 

Solution of the elastic problem. Let us denote by 
the stress tensor which is related to the strain field e,fc 
by the usual Hooke's law. Then, the stress in the new /3- 



phase can be written as a. 
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unchanged in the mother ev-phase: a 
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tensor is related to the lattice strain e° ik by Hooke's 
law 
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where E is the Young's modulus and v is the Poisson 
ratio. Thus, the mechanical equilibrium conditions at 
the interface require us to introduce the interface force 
density in the equilibrium equation d&i k /dx k = /j, 



fi = ip. 
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(5) 



where the normal vector n points from the /3-phase into 
the a-phase. Because the forces act only at the interface, 
the displacement field can be written as an integral over 
the interface surface: 

«i(r)= f G ik {v,v')f k {v')dS' , 



where Gik{v, r') is the so called Green's tensor (see, for 
example, [16]). Then, the strain field is determined by: 

If ( dG km (r,r') dG im (r,r') \ , , 
Qfe(r) = 2 J { d Xi + d Xk ) /m(r )dS - (6) 

This strain is fully defined by the corresponding Green's 
tensor and Eqs. (4,5). The strain components e rr , e ST , e ss 
are continuous across the interface while the other com- 
ponents are discontinuous: 
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Taking these jumps of the strain field at the interface into 
account, one can find the elastic contribution to the local 




FIG. 1: Steady state growth of a bicrystal (/3/f3'). The struc- 
ture propagates with a constant velocity v along the y axis. 
The /3-phase corresponds to 9 — 2ir/3 and the /3'-phase cor- 
responds to 9 = —2n/3. 



equilibrium condition, Eq. (3). 
forward calculation leads to: 
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where is the strain in the a phase at the interface. 
Note that the expression above is a complicated intcgro- 
differcntial functional of the interface shape. 

Diffusional growth. For simplicity, we consider transi- 
tions in pure materials and assume that the heat diffusion 
constants are equal in both phases (the so-called symmet- 
rical model). Without loss of generality, we assume that 
the (3 phase is the low temperature phase. We introduce 
the dimcnsionless temperature field w = c p (T — T OQ )/L, 
where L is the latent heat, c p is heat capacity, and is 
the temperature in the a phase far away from the inter- 
face. The temperature field w obeys the following heat 
diffusion equation and boundary conditions 

DV 2 w = dw/dt, (8) 
v n = Dn(Vwp\ int - Vw a \ int ), (9) 
w\ mt = A-d oK + T eq c p 5F el /L 2 7 (10) 

where do = jT eq c p / ' L 2 is the capillarity length, k is the 
curvature of the interface, which assumed to be positive 
for convex interfaces, D is the thermal diffusion constant, 
and T eq is the equilibrium temperature for the flat inter- 
face without elastic effects i.e, it is determined by the 
condition F°(T eq ) = F^(T eq ). We also introduce the 
dimcnsionless undercooling A = c p (T eq — T^/L. The 
physics underlying Eqs. (8-10) is quite simple. The in- 
terface moving with the normal velocity v n releases the 
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latent heat. Requirement of heat conservation at the 
interface gives Eq. (9). The local thermodynamical equi- 
librium at the interface, Eq. (3), implies Eq. (10), which 
gives the equilibrium value of the temperature at the in- 
terface taking into account the curvature corrections and 
the elastic effects. The thermal field can be eliminated 
by using the Green's function techniques (see for example 
[14]), and consequently together with a proper Green's 
tensor Gjfc(r, r') for the elastic field, the set of Eqs. (6- 
10) can be incorporated into a single integro-differential 
equation for the shape of the solid-solid interface. 

Dilatation. First, we consider the dilatational case, 
e ifc = £ $ik, where the bond lengths of the new phase are 
uniformly longer or shorter in all directions in comparison 
to the original phase. In this case, the elastic contribu- 
tion to the local equilibrium condition at the interface, 
SF( el ) = — e 2 E/(l — v), is a constant along the interface 
for any interface shape. This result corresponds to the 
elastic hysteresis mentioned above and it can be obtained 
using the analogy of this elastic problem to the problem 
of thermal expansion for a given temperature field [16]. 
As a consequence, the equilibrium interface temperature 
is shifted by a constant value, and the problem is equiv- 
alent to the problem of the classical dendritic growth. 
The anisotropy of the surface energy is a prerequisite for 
the existence of the dendritic solution in this case, as 
we have already mentioned above. Note that this exact 
result is valid only under the assumptions of our model 
i.e., the isotropic elasticity and equal elastic moduli in 
both phases. It also serves as a nontrivial check of our 
numerical code. 

Shear transition. Let us consider now a simple type 
of transition in hexagonal crystals involving the shear 
strain. For the transitions lowering the symmetry from 
Cq to C2 the shear strain in the basic plane appears. This 
is the case, for example, in hexagonal-orthorhombic tran- 
sitions in fcrroclastics (see [17] and references therein). 
Let the principal axis Cq be orientated in the z direc- 
tion. We assume that the system obeys a translational 
invariance in this direction, and thus it is effectively two- 
dimensional. By proper choice of the crystal orientation 
around the main axis in the initial phase, we obtain the 
new phase in three possible states, having the following 
nonvanishing components of the strain tensor e? fe : 

e° xx = -e yv = ecos26 e° xy = e sin 20, 



(11) 



where the angle 9 has three possible values: 9 = 
0, ±27r/3. The situation of 9 = corresponds to the 
single crystal growth and 9 = ±2ir/3 to the bicrystal 
growth (Fig.l). Because the elasticity of hexagonal crys- 
tals is axisymmetric in the harmonic approximation and 
(-i Z = £iz = 0, we can use the isotropic theory of elastic- 
ity i.e., expressions for free energy (1,2) remain valid [16]. 
The moduli of the effective isotropic elasticity, A and /z, 
can be expressed in terms of the elastic constants of the 
original hexagonal crystal. The two-dimensional Green's 



tensor for isotropic materials is [18] 

1 + V ( X{Xk 
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Eliminating the thermal field, we obtain the steady state 
equation for the shape of the solid-solid interface. In the 
co-moving frame of reference, this equation reads: 
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where rj = [(x — x') 2 + (y{x) — y(x')) 2 ]i , and K is the 
modified Bessel function of third kind in zeroth order, 
and p = vR/2D is the Peclet number. All lengths are 
reduced by the radius of the curvature R of asymptotic 
Ivantsov parabola. 

In the asymptotic region (\x\ — ► 00), the strain e-sfc 
decays and the local contribution to 5F el (second term in 
Eq. 7) has a constant value, 6F el = -E(e yy ) 2 /2(l - v 2 ). 
It follows from this relation that the temperature shift 
(clastic hysteresis) for the growth of a single crystal is 
four times larger than for a bicrystal. Although both 
bicrystal configurations [J3j (3' and j3' / (3) are energetically 
equivalent far from the tip, the symmetry is broken by 
the choice of the propagation direction. Therefore, in the 
following we will discuss the most favorable configuration 
of bicrystals [12], as presented in Fig. 1. 

Let us introduce the shifted, due to the elastic hystere- 
sis, undercooling: 

A = A — A e (, A el = T eq c p Ee 2 /%{l-v 2 )L 2 . 

The dimensionless parameter A e ; describes the strength 
of the clastic effects. The relation between this 
shifted undercooling A and the Peclet number is given 
by the two-dimensional Ivantsov formula [3]: A — 
^/pnexp (p) crfc( % /p). 

The presence of the twin {(3/13') boundary leads to 
additional effects. First, calculating the strain field, 
Cik, which enters in the expression for 5F^ el \ in ad- 
dition to the integrals along a/ (3 and a/ (3' interface, 
the integration should also be performed along the twin 
boundary /?//?'. The force density at this boundary is: 
fx = 0, f v = £'v / 3/(l + v). Second, the equilibrium an- 
gle 4> & t the triple junction (Fig. 1) is given by Young's 
law: 27 sin = 7b, where 7 is the surface energy of the 
(3 1 a. interface and 76 is the surface energy of the twin 
boundary. 

Eq. (13) is a complicated nonlinear integro-differential 
equation for the interface shape. We should find a solu- 
tion of this equation which has a proper angle </> at the 
triple junction and which is close to the Ivantsov parabola 
(y = —x 2 /2) in the tail region. Note that without elastic 
effects, this problem is equivalent to the classical den- 
dritic growth problem with isotropic surface tension. The 
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FIG. 2: Stability parameter a* versus A e i/p for two values of 
0: the dashed line corresponds to <j> — 0, and the solid line 
corresponds to = ir/6. 



latter does not have a solution with angles > [1, 2]. 
This statement can be expressed in the following form. 
For any given positive values of the Pcclet number p and 
the so-called stability parameter a = d /pR, the sym- 
metric solution which is close to the Ivantsov parabola in 
the tail region has an angle at the tip = f(a,p) < 0. 
The limit a = and = is a singular limit for that 
problem. For example, Meiron [19] calculated the angle 
as a function of a for several values of the Peclet num- 
ber with isotropic surface tension numerically and found 
that the angle < for any positive a. 

Now, we discuss the numerical results obtained by the 
solution of Eq. (13) in the spirit of Rcf. [19]. In the 
important regime of small Pcclet numbers, the eigenvalue 
a = a* (0* , A e i , p) depends only on the ratio A e i jp for a 
fixed angle 0. While the strength of the elastic effects 
is assumed to be small, A e ; <C 1, the control parameter 
A e i lp can be varied in a wide region in the limit of small 
p. The eigenvalue a* as a function of A e i jp for two values 
of the angle, = and = tt/6, is shown in Fig. 2. The 
situation with ss is realized if 7b -C 7, while ks tt/6 
corresponds to 76 ~ 7. The Poisson ratio was fixed to 
v = 1/3. 

The most remarkable feature of these results is that 
we do find dendritic solutions for the isotropic surface 
tension in the presence of the elastic effects. In this sense, 
the elastic effects serve as a new selection mechanism. 
We note that a* becomes large for large values of A e i /p, 
while in the classical dendritic growth a* is always small, 
being controlled by tiny anisotropy effects. Thus, the 
growth velocity, v = 2Da*p 2 /do, can be much larger 
due to elastic effects, compared to the classical dendritic 
growth. 

In the case = tt/6, the solution exists only beyond 
some critical value of the control parameter, A e i/p. The 
lower branch of this solution has a negative curvature 



near the triple junction and it is presumably unstable. 
The general structure of the theory suggests that in the 
case = 0, the curve should start from the origin (a = 
0,A e i/p — 0) which is a singular point of the problem. 
However, numerics becomes very difficult in the vicinity 
of this point. 

We have also performed several runs for the single crys- 
tal growth (0 = in Eq. (11)). As in the classical den- 
dritic growth, we have not found solutions with a smooth 
tip, = 0. Moreover, for negative values of 0, where so- 
lutions exist and correspond to the growth along a grain 
boundary in the mother a phase [20] , the selected stabil- 
ity parameter a* is a decreasing function of A e i/p. This 
is in strong contrast with the results for the bicrystal 
growth. 

We hope that our results will stimulate new experimen- 
tal and theoretical investigations in this interesting field. 
On the theoretical side, a challenging problem would be 
the investigation of the doublon patterns in the case of 
bicrystal growth. With this scenario, two crystals of the 
f3 phase would be separated by a thin film of the mother 
a phase and not by the twin boundary. This pattern 
is a competitive structure at least in the range of the 
small control parameter A e i/p, where the dendritic pat- 
tern with the twin boundary does not exist. 
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